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We point out a simplicity that arises when we use an interaction in which only an energy with 
odd J is non-zero. The emphasis is on J = Jmax and in particular J = 9"'' in the 39/2 shell. It 
is noted that high overlaps can be deceptive. In many cases a single set of unitary 9-j coefficients 
gives either an exact or a very good approximation to the wave function of a non-degenerate state. 
The many degeneracies that occur in these calculations are discussed and explained. 
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I. INTRODUCTION 



The purpose of this work is to study the properties of 
a very simple interaction in a single-j-shell model space 
of neutrons and protons. A single proton-neutron (pn) 
pair in this space can have a total angular momentum 
from J = to J = Jmax = 2j. The even- J states have 
isospin T = 1, i.e. they are members of an isotriplet — 
there are analog states of two neutrons and of two protons 
with these even angular momenta. The odd- J states have 
isospin T = 0, i.e. isosinglet — in this model space, they 
only are present in the pn system. 

There has been much study in journals and in text- 
books of the J = 0, T = 1 "pairing interacton", i.e. where 
only J = two-body matrix elements are non-zero (and 
attractive) [H-O- Although the pairing interaction is not 
realistic, the consequences of these studies has yielded re- 
sults whose importance was well beyond expectation. Ex- 
amples are seniority classifications, reduced isospin and 
indeed these studies can be regarded as precusors to the 
BCS theory in condensed matter theory. In turn the BCS 
theory is important in the nuclear context for explaining 
moments of inertia of deformed states. In the pn system, 
it is known that not only the J = two-body matrix 
element lies low but also J = 1 and J = Jmax — 2j. As a 
counterpoint to the J = paring interaction, we will here 
consider J = Jmax pairing interaction, for which V{2j) 
is non-zero and attractive and all other two-body matrix 
elements are zero. Such an interaction only acts between 
a neutron and a proton, not between two identical parti- 
cles. This study will reveal some surprising results. 

It should be noted that there has been a recent flurry of 
interest in Jodd pairing and in particular the case where 
Jodd is equal to Jmax- In a receirt work, C. Qi et al. [S] 
proposed a wave function for ^®Cd in which each of two 
pn pairs couple to maximum angular momentum, which 
for the (7g/2 shell is J = 9. More generally, they con- 
sider the overlap of a wave function in which one pn pair 
couples to Ji and the other to J2, with basis states of 
two protons coupling to Jp and two neutrons to J„. In 
the latter basis, the effect of the Pauli principle is much 
clearer. Once one has antisymmetric wave functions of 
two protons, which is achieved by limiting the angular 
momenta to even Jp, and likewise the two neutrons to 



even J„ , one has a wave function which satisfies the Pauli 
principle. For the former wave function [(pn) Ji {pn)J2Y, 
one does have to antisymmetrize and a priori thiirgs look 
complicated. The authors, however, make the statement 
"The overlap matrix automatically takes into account the 
Pauli principle". Here the overlap matrix element is 

0^{JiJ2;JpJn) = {[pn{Ji)pn{J2)Y\[pp{Jp)nn{Jn)Y) , 

(1) 

that is, the unitary 9-j symbol, which we will here call 
U9-J. 

We here note that work we previously did shows that 
certain U9-j coefficients form components of the J = 0+ 
ground state wave function of a Jodd pairing interaction. 
We here expand on this work. We feel our method is 
very simple to understand. We refer to a previous work 
by E. Moya de Guerra et al. [T] and explicitly to Eqs. 
(70) and (74). In that work we describe the wave func- 
tions in a [{pp)Jp (nn) J„] basis. The Pauli principle is 
easily satisfied by constraining Jp and J„ to be even. 
Our previous example was ^"^Ti, but the same mathe- 
matics holds for ^^Cd. We previously considered vari- 
ous schematic interactions as well and the more realistic 
MBZ interaction taken from experiment [8], for which de- 
tailed wave funtions were subsequently published in the 
archives [21 (with some modification of the two-body ma- 
trix elements). In this work we use the symbol I for the 
total angular momentum of the state and J otherwise. 



II. CALCULATIONS 

We now proceed to the calculations. We will be exten- 
sively using these two relations for the [/9-j's: 
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where s = J24 + J23 - J34 - 1. 

Let us first consider the 7 = 0+ states in ^^Cd (or 
^*Ti). For most interactions, the diagonahzation is a 
fairly complicated procedure. However for certain inter- 
actions it is much easier. For example, the interaction 
used in Ref. 01 was one in which all two-body matrix el- 
ements were set equal to zero except for the J = 1, T = 
two-body matrix element. In such a case, the matrix el- 
ement of the secular four-particle Hamiltonian factorizes 
(the same result holds for any odd- J interaction). This is 
the key point. The basis states are [Jp JpY^°. We have 

Hj^,j^, =E{JoM)f{Jp)f{Jp'), (4) 
where f{Jp) is twice the U9-j symbol: 

{ j j Jp] 

= 2{2Jp + l)(2Jodd + 1) < J J Jp > (5) 

[Jodd Jodd J 

If we write the wave function as ^ Xj^j [JpJpY^^ (as 
in Ref. Q), then it was shown in Ref. Q that Xj^j^ 
is proportional to f{Jp). The other eigenstates are de- 
generate and, if E{Jodd) is negative, they are at higher 
energies. In other words, what we have shown in Ref. [7| 
is that the wave-function components Xj^j^ of the low- 
est / = 0+ state are proportional to the overlap factor of 
Ref. [gI; alternately, they are equal within a normaliza- 
tion to the J79-j coefficients. 

The eigenvalue is given by 

E{I = 0+) = i?( Jodd)| ^ f{Jp)Xj^j^ ? (6) 

Note that our very simple interactions are charge inde- 
pendent. This means that the lowest (non-degenerate) 
7 = 0+ state has good isospin, presumably T = 0. It is 
amusing that we can assign the isospin quantum number 
to a wave function with C/9-j coefficients. 

In Table |T] we present the wave functions for the fol- 
lowing interactions: 

CCGI: Ref. [lo!|. 

E(0): V{0) = -2.0000 MeV; aU other matrix elements 
are zero. 

E(9): V{0) = -2.0000 MeV; aU other matrix elements 
are zero. 

E(0,9): y(9) = £^(0) = -2.0000 MeV; all other matrix 
elements are zero. 

E(l): V'(l) = -2.0000 MeV; all other matrix elements 
are zero. 

We consider the second column (CCGI) as the real- 
istic interaction to which the other interactions should 
be compared. The simplest thing we can do is give the 
overlaps of the above interactions with CCGI. They are 



Table I: Wave Functions of the 7 = ground state of Cd 



for various interactions. 




CCGI 


E(0) 


E(9) 


E(0,9) 


E(l) 


^00 


0.7725 


0.8563 


0.6164 


0.8103 


0.2903 




0.5280 


0.1741 


0.7518 


0.4814 


0.5704 


X44 


0.2915 


0.2335 


0.2385 


0.2514 


0.5190 




0.1704 


0.2807 


0.0233 


0.1718 


0.1586 


^88 


0.1020 


0.3210 


0.0005 


0.1831 


-0.5540 



respectively 0.9020, 0.9467, 0.9944, and 0.6484. We find 
that E(9) gives higher overlap than the much studied 
E(0) pairing interaction and a much higher overlap than 
E(l) . This might lead one to believe that the idea of 
J = 9+ pairing is a valid concept. But overlaps can be 
deceiving. We also present E(0,9), where the only non- 
vanishing matrix elements are for J = 9+ and 0+, both 
set to —2.0000 MeV. Now the overlap is even higher — 
0.9944. This might not be startlingly different than 
0.9467, but let us now look at the energies of the low- 
est even-7 states in Table [TTl They are given respectively 
for interactions CCGI, £^(9)), and E(0,9). The results 
of the second column (CCGI) were previously given [ll| 
and the point was made that the 7 = IG"*" state is iso- 
meric since it lies below the lowest 14+ and IS"*" states. 
This is in agreement with experiment (l2| . 



Table II: Calculated spectra of yrast even-7 states in Cd for 
above mentioned interactions. 



r 


CCGI 


S(9) 


E(0,9) 


0+ 


0.0000 


1.0587 


0.0000 


2+ 


1.0812 


1.0589 


1.2740 


4+ 


2.1096 


1.0588 


1.8584 


6+ 


2.8883 


1.0588 


2.3929 


8+ 


3.2302 


1.0571 


2.5125 


10+ 


4.8815 


1.0464 


3.2142 


12+ 


5.3394 


0.9670 


3.1348 


14+ 


5.4031 


0.6570 


2.8247 


16+ 


5.2247 


0.0000 


2.1678 



We see that despite the 0.9467 overlap, the even-7 spec- 
trum for Intl in which only the J = 9+ matrix element 
is non-zero is drastically different than CCGI. First of 
all, the ground state does not have 7 = 0, rather it has 
I ~ Jmax = 16 and indeed the two spectra seem to have 
nothing to do with each other. 

Let us briefly digress and look at the spectrum for Intl 
for its own sake. It is quite remarkable. The energies of 
the 7 = 0,2, 4, 6, and 8 states are very close to each other, 
differing at most by 0.002 MeV and the 7 = 10+ state 
is 0.012 lower. All six states are essentially degenerate. 
Then there is a drop in energy with 7 = 16+ becoming 
the ground state. Such a strange spectrum and this for 
an interaction that gives a 0.9467 overlap with a realistic 
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interaction for the 1 = 0"'" state. 

In the last column of Table |TT1 we improve things by 
also lowering the J = 0+ matrix element to the same 
value as for J = 9+, -2.0000 MeV. The spectrum is 
better, with 7 = 0+ now the lowest state, but it is far 
from satisfactory. Even an overlap exceeding 0.99 does 
not guarantee overall good results. Clearly all two-body 
matrix elements come into play. 

As noted above, the eigenfunction of the low- 
est 1 = 0+ state for the £'(9) interaction is 
N{{fY''{fy"\{ff{ffy- It can be shown that the 
normalization factor is \/2. 

For the 7=1+ states with the i?(9) interaction, the 
secular matrix is also separable. This is not true for other 
values of Jodd • If we were to replace / = by / = 1 in the 
above expression, all the t/9-j coefficients would vanish. 
We must make a different choice. The eigenfunction of 
the lowest 7=1+ state is then given by a single set 
of C/9-J coefiicients: 2{{j'^Y'-{fY-\{pf{j'^fY=^. This 
state has isospin T = 1. Indeed all / = 1+ states in this 
model space have isospin T = 1. The other four 7=1+ 
states are degenerate at a higher energy. 

For states with / = 2 or higher, the secular matrix 
is no longer separable — rather it is a sum of separable 
terms. The eigenvalue equation is 

^Y.^{fYYfY-m\jy^Y^ 
X E m'^'{fY-'m'{3y^YD{j,..j,,) = 

- XD{Jp,Jn) (7) 

For 7 = 2+ there are two terms corresponding to J^; ~ 7 
and 9; for 7 = 3+ the values are ~ 6 and 8, etc. 

Despite the complexity of the above equation, 
there are some surprising results. The eigenfunc- 
tion components of the lowest 2+ state are numer- 
ically extraordinarily close to the single U9-j sym- 
bols V2{{YYYfy^"\{YY{YYY"'^- Furthermore, the 
next 2+ state has also components exceedingly close to 
'^HYy-^Yf y^^liYYiYY)^^^- This is by no means obvi- 
ous because, as mentioned above, the interaction involves 
a sum of two separable terms corresponding to J^, = 7 
and 9. We can explain this result by performing the sum 
over even Jp and even J„. We first note schematically 

even JpJn 

= E + E(-i)'^ + E(-i)'" + E(-i)''^'" (8) 

The first term vanishes because of Eq. ^ . In the last 
term one of the J79-j's has two rows that are the same, 
which means that the only non-vanishing terms in the 
sum have (Jp + Jn) even. Thus, the last term is the same 
as the first term — zero. The two middle terms are the 



same, so we get 

E - ^E(-i)'''((j'')'^(j'')'"K^'')'i(j'')')'='x 

even JpJn 

^{{.fYYfY-\{fY{.fyY='- 
= -liifYifYKfYifYY-'- (9) 

We obtain the above by using two orthogonality relations 
for 9j-symbols as shown e.g. in Nuclear Shell Theory 
(p. 516) by de-Shalit and Talmi ||4|]. We call the right- 
hand side of Eq. ^ the overlap. 

Using similar arguments, one can show that the nor- 
malization for the [9,9], N, is such that 

= l-\{{fnjr\{fYifYY-'= (10) 

= I- ^0.00001209813 = 0.499993950935 
2 2 

For the [9 , 7] case, we obtain 

^ \-\{{fnjy\{fY{fYY-'= (11) 

= T + ^0.00075253477 = 0.250376267385 

4 2 

To get this latter result, we use the following relation- 
ship 

(12) 

Therefore, we obtain that the overlap is exceedingly 
small for the gg/2 shell. From Eqs. pUl) and PT|) . we find 
that the normalizations are 1.414222 and 1.998497, the 
latter slightly smaller than 2. The overlap in Eq. ([9]) is 
0.00009113 and, if we include the normalization factors, 
we get 0.00025756. 

In lower shells the deviations are larger. For exam- 
ple, in the ^5/2 shell, we replace [9,9] by [5,5] and 
[9,7] by [5,3]. The overlap is now 0.0107, small but 
not zero. The norms are no longer Y2 and 2, but rather 
1.4161 and 1.9204. With normalized states the overlap 
is 0.0290. The corresponding numbers in the f'j/2 shell 
are 0.001037, 1.41434, 1.9875 and 0.002916. The overlap 
is here smaller than in the shell, but it is not zero. 
As one can see above in the 39/2 shell, the overlap is an 
order of magnitute smaller than in the /y /2 shell. We can 
speculate that the overlap might vanish in the infinite-j 
limit. 

We can see in Table IIIII that the results for matrix di- 
agonalization for both 7 = 2+ states yield wave function 
components which are very close to the normalized U9-j 
coefficients. In fact, they are so close that one could won- 
der if they are exactly the same. But they are not. As 
seen in Eq. (|9]) , the two U9-j sets corresponding to [9 , 9] 
and [9 , 7] are very nearly orthogonal, but not quite. 

The above 7 = 2 states have isospin T = 0. The E{9) 
interaction also yields a T = 1 non-degenerate state with 
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components '2,{{ijY{jjY'\{jjy^''{jjY'^)'^ . This is a pure 
state — it does not mix with any other T — 1 state. This 
is because there is only one way of forming an / — 2, 
T = 1 state from UQ-j symbols, i.e only one possible J^- 

For / = 3, T = there is also a pure state 
'^{{ji)'^{ijy\{jiY^{iiY")^- This wave function changes 
sign under the interchange of Jp and J„ . It cannot admix 
with a state with Jx = i for which there is no change of 
sign when Jp and Jn are interchanged. 



Table III: Comparison for the first two 7 = 2+ states of the 
matrix diagonalization with the i?(9) interaction and with 
normahzed U9-j components. We give the energy in MeV in 
the second row. 



Jj,,J„] £(9) m-j E{9) U9-j 

1.069 3.0558 

[0,2] 0.5334 0.5338 0.1349 0.1351 

[2,2] -0.4707 -0.4708 0.5569 0.5567 

[2,4] 0.3035 0.3035 0.3188 0.3189 

[4,4] -0.1388 -0.1390 0.6300 0.6299 

[4,6] 0.0531 0.0531 0.1320 0.1320 

[6,6] -0.0137 -0.0138 0.1350 0.1350 

[6,8] 0.0025 0.0025 0.0114 0.0114 

[8,8] -0.0003 -0.0003 0.0052 0.0052 



It turns out that all the other lowest even- J states 
have eigenfunctions close although not exactly equal to 
V2{{fy''ify^"\if fif)y. in Table HV] we compare, 
as an example, the wave function of the J = 8+ state. 
In the second column, we give the single U9-j symbols 
(normalized) and in the third column we give results of 
diagonalizing the E{9) interaction. Since the coefficient 
[Jp, Jn] is the same as [J„, Jp], we list only one of them. 
The overlap of the two wave funtions, (tAi, V'a), is 0.9944. 



Table IV: Comparing the wave functions of a single U9-j sym- 



bol with Jx ~ 9 with a full diagonalization of -E(9) for the 
lowest 1 = 8+ state in '*'^Cd. 

[Jp , Jn] lJ±j E{9) 

[0 , 8] 0.0630 0.0644 

[2 , 6] 0.4299 0.4271 

[2,8] -0.0522 -0.0513 

[4 , 4] 0.7444 0.7456 

[4,6] -0.1803 -0.1729 

[4 , 8] 0.0256 0.0280 

[6 , 6] 0.0521 0.0657 

[6,8] -0.0076 -0.0012 

[8,8] 0.0011 0.0047 



III. DEGENERACIES 

With this E{9) interaction, we get several degenerate 
states with an absolute energy zero. In some detail, for 
7 = there are five states, three with isospin T = and 
two with T ~ 2. There is one non-degenerate state at an 
energy 2V{9) {V{9) is negative). The other four 7 = 
states have zero energy. For 7 = 1 all states have isospin 
T = 1. There is a single non-degenerate state at ^(9), 
the other three have zero energy. For 7 = 2 there are 
twelve states — six have T = 0, four have T — 1 and two 
have T = 2. There are two non-degenerate T — states 
with approximate energies 21^(9) and V{9) respectively, 
and one non-degenerate T = 1 state with energy V{9). 
The other nine states have zero energy. To understand 
this, take a wave function 

and the corresponding energies 7?" = (^'"77^"). Con- 
sider the sum V E". We have 

^C"(Jp, J„)C"(Jp,, J„0 = <5j,.j^,5j„.j„, (13) 

a 

Thus 

- H {[JpJn]'H[JpJnY)= (14) 

a JpJn 



= 4n9)EEKj^')''''(j^')''"i0'j')'0j) 




JpJn J A 

even 

This expression does not depend on the detailed wave 
functions. Using the properties of t/9-j's, we can show 
that Ea-S" 2F(9) for 7 = 0, V{9) for 7 = 1, and 
4V^(9) for 7 = 2. But we can alternately show, using the 
explicit wave functions, that for 7 = the energy of the 
lowest state is 2^(9). Hence, all the other states must 
have zero energy. A similar story for 7 = 1. The 7 = 2 
state is a bit more complicated because of the coupling 
between two states, however small it is. Still one can work 
it through and see that the 4F(0) energy is exhausted by 
the two T = and the one T = 1 non-degenerate states. 

IV. CLOSING REMARKS 

In closing, we note that the subject of J,nax pairing 
is currently a very active field. Besides the work of Qi 
et al. [6|], t here are related works by Zerguine and Van 
Isackerjl^. Cederwall et al. |il4] and Xu et al. [15]. The 
topic of J-pairing interactions has also been addressed 
by Zhao and Arima [16]. In this work we expand on our 
2003 work Q by renoting that the Hamiltonian matrix 
for a 2p-2n system for 7 = 0+ states in a single j-shell 
is separable for a simple interaction which is non-zero 
only for a single odd angular momentum. This leads 
to an eigenfunction with components proportional to a 
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single set of unitary 9-j symbols. We apply this to the 
J = Jmax interaction. The single set of U9-j components 
form the eigenfunction not only for the lowest 7 = 0+ 
state, but also of the lowest 7 = 1+ state and, to a 
surprisingly excellent approximation, for the lowest two 
7 = 2+ states. A single set of U%-j coefficients yields a 
good approximation for the higher yrast even-/ states. 
We use the J = Jmax interaction to confirm the obser- 
vation of Ref. that the resulting wave function has 
a fairly high overlap with that of a realisitc interaction. 



although we note the energies are not realistic. We also 
note that this wave function for the lowest 7 = 0"'", l"*", 
and 2+ states has good isospin. We have other examples 
too: we have found a quantum number (see Eq. ([7])) 
which can, either exactly or approximately, help classify 
some of the states. 

Our cautionary remarks are for the topic of pairing in 
general, but are not intended specifically for any of the 
works mentioned in this paper. 
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